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I.  INTRODUCTION 

The  most  important  aspect  of  pump  design  is  the  design  of 
the  Impeller.  It  is  this  component  that  will  determine  how  well 
the  pump  performs,  whether  or  not  it  will  meet  the  conditions 
required  of  it.  In  the  design  problem  the  desired  flow  condition 
Is  specified,  and  it  is  required  that  the  impeller  blading  provide 
the  proper  fluid  turning  without  absorbing  excessive  power  due  to 
fluid  losses  within  the  blade  row.  It  is  therefore  imperative 
that  the  impeller  be  designed  as  accurately  as  possible,  incor- 
porating into  the  design  sound  fundamental  hydraulic  principles. 

The  three  basic  types  of  impellers  — radial  flow,  axial  flow, 
and  mixed  flow,  are  illustrated  in  Figure  1.  The  blades  of  purely 
radial  flow  impellers  are  designed  by  way  of  tangential  arcs. 

In  this  method,  illustrated  in  Figure  2,  a blade  section  is 
described  by  circular  segments  which  smoothly  increase  the 
tangential  or  circumferential  sweep  of  the  blade  as  the  radius 
is  increased  from  inlet  to  exit.  The  blades  of  axial  flow 
impellers  are  designed  by  considering  cylindrical  sections 
at  a constant  radius.  Here,  blade  shapes  are  laid  out  for  each 
radial  location  selected,  in  much  the  same  manner  as  a cascade, 
increasing  the  angular  wrap  of  the  blade  as  the  axial  distance 
is  traversed  from  inlet  to  exit.  Figure  2.  In  each  of  these  cases 
the  designer  is  required  to  work  only  with  two  dimensions,  the 
other  being  held  constant  or  arbitrarily  specified  for  each  blade 
section. 


The  mixed  flow  Impeller  as  well  as  the  Francis-type  impeller, 
both  of  which  have  predominantly  axial  inlet  and  radial  exit  flow. 
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do  not  lend  themselves  to  either  of  these  two  methods.  They  are 
truly  three-dimensional  problems  in  that  the  blade  shapes  must 
change  in  the  axial,  radial,  and  tangential  direction. 
Unconventional  axial  flow  impellers,  where  the  hub  radius 
Increases  from  inlet  to  exit,  are  also  three-dimensional  in 
nature  since  the  streamlines  shift  outward  as  they  travel  through 
the  blades.  A method  to  handle  these  complicated  three-dimensional 
design  problems  is  presented  in  this  report. 

All  pump  designs  begin  by  utilizing  simple  one-dimensional 
considerations,  which  specify  the  required  parameters  at  the 
pump  inlet  and  exit  stations.  More  sophisticated  designs  employ 
two-dimensional  considerations  to  analyze  the  axlsymmetric  flow 
field  in  an  attempt  to  determine  what  takes  place  between  the 
pump  inlet  and  exit.  However,  the  design  of  many  pump  impellers, 
as  has  been  pointed  out,  is  truly  a three-dimensional  problem, 
and  ways  must  be  developed  to  consider  the  real  fluid  flow. 

The  aim  of  each  of  these  efforts  is  to  somehow  construct  impeller 
blades  which  will  allow  the  pump  to  operate  as  designed. 

This  report  is  meant  to  guide  the  designer  through  the 
design  process  with  the  help  of  the  computer.  The  designer  will 
first  consider  a one-dimensional  approach  to  pump  design  and 
attempt  to  optimize  the  pump  parameters  to  satisfy  performance 
and  operating  requirements.  Next,  the  proper  two-dimensional 
axlsymmetric  pump  geometry  is  selected  for  this  particular  design. 
Then,  in  the  most  Important  step,  the  blade  shapes  in  three- 
dimensional  space  are  designed  and  constructed. 


F 
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This  process  If  undertaken  with  the  assistance  of  a computer 
program  written  in  the  Fortran  computer  language.  The  computer 
Is  operated  In  conjunction  with  a CALCOMP  plotter  which  constructs 
all  of  the  drawings  that  are  normally  done  by  hand  at  the  drafting 
table.  Because  of  the  speed  of  the  digital  computer,  a design 
can  be  completed,  modified,  and  redrawn  several  times  as 
required;  a task  that  would  take  the  engineer  and  draftsman  weeks 
to  accomplish. 

The  method  of  design  used  by  the  computer  is  that  of  George 
F.  .Wisllcenus  as  outlined  in  Reference  9.  It  makes  use  of  a 
conformal  mapping  which  lays  out  the  streamsurfaces  on  a 
rectangular  grid.  In  this  way,  two-dimensional  blade  shapes 
can  be  constructed  which,  although  distorted  in  the  mapping,  are 
true  representations  of  the  blading  in  three-dimensional  space. 

The  output  of  this  method  provides  for  the  mechanical  construction 
of  the  blading  by  producing  blade  shapes  described  by  board 
sections  which  can  be  built  up  layer  by  layer,  spanning  the 
entire  impeller  channel. 

By  using  the  computer,  steps  which  would  be  unnecessarily 
time-consuming  by  hand  can  be  undertaken.  In  this  way,  many 
more  data  points  can  be  used  to  ensure  a more  accurate  design. 

The  computer  program  was  directly  adapted  from  the  graphical 
procedure,  and  several  refinements  in  the  original  design  method 
were  incorporated  in  its  development.  The  program  is  capable  of 
employing  techniques  characteristic  of  the  Mean  Streamline 
Method  of  blade  design^  which  has  been  used  with  much  success 
in  the  design  of  conventional  axial  flow  machines. 

I 

I 

i 

i 
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An  existing  design  of  a Francis-type  Impeller  Is  Included 
as  an  example.  This  Impeller  represents  one-half  of  a double- 
auction  centrifugal  pump,  five  of  which  would  be  located 
transversely  on  a high-speed  surface  planing  craft.  This 
configuration  makes  use  of  a wide  aspect-ratio  Inlet  to  Ingest  a 
large  portion  of  the  boundary  layer  flow,  thereby  reducing  the 
amount  of  internal  diffusion  required  to  slow  the  fluid  velocity 
to  that  of  the  pump  inlet  velocity.  The  advantages  of  this 
system  are  a reduction  in  the  wetted  system  vreight  and  a decrease 
In 'the  momentum  of  the  incoming  flow. 
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II.  OVERVIEW 

Since  the  entire  design  process  has  been  Incorporated  Into 
the  computer  program,  a discussion  of  the  design  method  can 
then  be  presented  by  examining  the  program  In  detail.  The  com- 
puter program  consists  of  seven  separate  modules  which  are  used 
consecutively.  The  designer  In  this  way  can  stop  at  any  point 
In  the  design  process  and  go  back  and  make  changes  In  the 
design.  In  fact,  the  first  two  modules  are  Intended  to  be  run 
several  times  In  order  to  optimize  the  pump  geometry  and 
operating  conditions.  Once  these  Items  are  fixed,  the  detailed 
three-dimensional  design  begins.  Each  module  plots  Its  results 
Independently  and  provides  data  to  be  used  as  Input  to  the 
remaining  modules. 

The  modules  used  In  the  computer  program  are: 

(1)  DESMOD  (Design  Module)  - which  allows  optimization  of  pump 
operating  and  performance  parameters,  sizes  the  Impeller, 
and  determines  the  velocities  at  the  Inlet  and  exit 
stations. 

(2)  GEOMTRY  (Geometry  Module)  - which  fixes  the  Inlet  and 
exit  stations  In  space  by  selecting  an  optimum  geometry 
for  the  axlsymmetrlc  Impeller  channel. 

(3)  SCM  (Streamline  Curvature  Method  Module)  - which  determines 
the  velocities  and  pressures  at  each  station  through  the 
pump  for  each  streamline  In  the  axlsymmetrlc  flow  field. 

(4)  STREAM  (Streamsurface  Module)  - which  divides  the  three- 


-*^1 


dimensional  streamsurfaces  Into  shapes  which  approximate 
squares  for  use  In  the  conformal  mapping. 
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BLADE  (Blade  Shape  Module)  - the  most  important  of  the 
seven,  which  constructs  the  actual  blade  shapes  on  the 
conformal  grid  (S  - 0 mapping). 

RADIAL  (Radial  Section  Module)  - which  transfers  the 
blade  shapes  to  an  R - Z mapping  on  the  profile  view  of 
the  Impeller. 

BOARD  (Board  Section  Module)  - which  plots  the  axial  or 
plan  view  of  the  impeller  illustrating  the  blading 
board  sections  on  an  R - 6 mapping. 


(5) 


(6) 

r 

i 

(7) 
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III.  DESIGN  MODULE  (DESMOD) 

The  one-dimensional  approach  to  pump  design  begins  with 
the  selection  of  the  design  operating  condition.  Required  as 
Input  to  the  design  are  certain  parameters  which  pertain  to  the 
overall  system  of  which  the  pump  is  a part.  These  system 
parameters  are:  the  head  rise,  H;  the  flowrate,  Q;  and  the  net 
positive  suction  head,  NPSH. 

The  type  of  impellers  considered  in  this  report  are  character- 
ized by  a change  in  tip  diameter  from  inlet  to  exit  station  and 
thus  are  a degree  more  complicated  than  a conventional  axial 
Impeller  design.  Therefore  in  order  to  help  in  understanding  the 
design  process,  a series  of  parameters  will  be  presented  which 
apply  to  any  type  of  pump®.  These  are: 
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where  the  hub-to-tlp  diameter  ratio  X is  Inputcd.  Since  there 
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are  more  unknowns  (seven)  than  equations  (five),  two  additional 
Inputs  are  required  at  this  point.  Typically,  the  flow 
coefficient  ((>  will  be  varied,  and  a sixth  equation  will  be 
added  to  determine  the  nondlmenslonal  net  positive  suction 
head  T as  a function  of  <p,  based  on  some  cavitation  criteria. 
Another  common  selection  Is  to  fix  the  shaft  speed  of  the 
pump,  normally  at  a convenient  speed  based  on  the  prime  mover 
and  its  reduction  gearbox.  In  either  case,  the  inlet  station 
is  designed  first,  based  on  a cavitation  requirement  and 
perhaps  optimized  on  the  flow  coefficient. 

These  equations  are  all  that  Is  needed  for  the  design  of 
an  axial  flow  pump  impeller;  the  subscripts  designating  the 
inlet  and  exit  stations  are  disregarded  In  this  case.  However, 
for  a mixed  flow  or  centrifugal  pump,  two  stations  of  different 
diameter  must  be  accounted  for.  The  normal  procedure  Is  to 
design  the  inlet  station  to  satisfy  a cavitation  or  system 
constraint,  and  then  to  design  the  exit  station  to  produce  the 
head  required  of  the  pump.  The  first  four  equations  are  used 
in  a manner  similar  to  the  axial  pump  design;  but  in  the  fifth 
equation,  the  head  coefficient  is  calculated  based  on  a tip 
blade  velocity  at  the  exit  station.  Thus,  a method  to 
determine  the  exit  station  diameter  or  the  diameter  ratio 

(D„  /D_  ) is  needed. 

2 ^1 

A direct  procedure  outlined  in  Reference  8 allows  one 
to  calculate  the  diameter  ratio  if  a value  is  selected  for 
the  retardation  of  the  tangential  component  of  the  relative 
velocity  (relative  to  the  impeller).  The  equations  are  as 
follows: 


n 
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where  an  estimated  hydraulic  efficiency  ri  and  an  expected  value 

n 


for  the  percent  prerotation  Xp  are  estimated.  These  equations 
can  be  combined  to  arrive  at  an  expression  for  the  diameter 
ratio: 
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where  the  diameter  ratio  has  been  substituted  for  the  tip 
velocity  ratio,  since  from  Equation  (A), 


U„  U„ 


D D 
U 2 


(] 


The  exit  station  is  also  characterized  by  its  dimension  b2. 


Figure  1;  thus,  another  input  is  required  to  complete  the  one- 
dimensional  design.  The  equations  needed  and  variables  are  as 
follows : 


Q (1  + Xj ) 
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and 


(14) 


where  Xy  is  inserted  as  the  reduction  in  the  meridional  component 
of  the  velocity  through  the  impeller,  and  is  the  estimated 
flowrate  leakage. 

At  this  point,  it  is  possible  to  construct  the  velocity 
diagrams  for  the  inlet  and  exit  stations  and  calculate  the  blade 
angles  required  of  the  impeller.  The  various  velocity  components 

and  fluid  angles  are  calculated  based  on  trigonometric  relations. 

* 

It  is  then  assumed  that  the  blade  angle  at  inlet  will  correspond 
to  the  inlet  flow  angle  6^^.  At  the  exit  station,  however,  the 
blade  angle  and  fluid  angle  will  not  be  the  same  since  the  real 
fluid  displays  an  irrotational  effect  resulting  in  what  can  be 
called  fluid  slip.  The  amount  of  slip  is  estimated  and  inputed 
as 


where  a star  (*)  Indicates  a quantity  for  which  the  fluid  angle 
and  the  blade  angle  are  assumed  to  be  identical.  The  exit 
station  velocity  diagram  is  then  constructed  for  both  starred  and 
unstarred  quantities.  A sample  output  from  DESMOD  is  listed  in 
Appendix  A.  The  inputs  and  outputs  are  shown  for  the  design 
condition,  and  velocity  diagrams  are  Included  to  help  visualize 
the  changes  in  the  flow  from  the  inlet  to  the  exit  station. 
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selected,  the  program  will  completely  construct  and  describe 
the  pump  geometry  without  any  further  Information.  For  case 
3,  It  is  possible  to  employ  the  required  geometrical  parameters 
In  such  a way  that  nearly  any  geometry  desired  by  the  designer 
can  be  specified.  For  the  designer  who  wants  to  investigate 
several  different  geometries.  It  Is  suggested  that  cases  1 and 
2 be  run  first.  Then,  If  necessary,  he  can  specify  the  geometry 
in  any  other  manner  by  running  case  3 as  often  as  desired. 

Each  streamline  is  constructed  by  using  two  circular  arc 
sections,  or  possibly  three  for  the  center  streamline.  The 
blade  leading  edge  is  also  defined  by  a circular  arc.  Thus, 
the  space  occupied  by  the  blading  is  identified,  the  trailing 
edge  being  defined  as  the  exit  station  at  the  tip  diameter. 

(For  case  3,  variable  Input,  a chamfered  or  oblique  trailing 
edge  may  be  used.)  Fictional  entrance  and  exit  stations  are 
located  1/2 (r_,  - r„  ) upstream  and  1/2 (b„)  downstream  of  the 

Tf  Hf  2 

actual  pump  configuration,  respectively.  These  sections  allow 
for  fluid  entering  the  leaving  the  axisymmetric  flow  field 
without  the  effects  of  curvature.  The  inner  and  outer  shrouds 
are  assumed  to  be  axial  at  inlet  and  radial  at  exit.  If  the 
two  circular  arcs  which  define  the  streamline  do  not  extend 
all  the  way  to  either  the  entrance  or  exit  section,  the  program 
will  do  so.  The  center  streamline,  however,  must  be  defined 
all  the  way  to  the  fictional  entrance  station.  This  is  the  reason 
for  allowing  it  to  consist  of  three  segments,  the  last  of  which 
may  be  a third  circular  arc  or  a linear  segment. 


j 

i 


The  data  are  generated  for  typically  30  stations  overall. 
Figure  A.  The  fictional  entrance  station  is  selected  as  number 
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8,  the  leading  edge  as  number  16,  the  trailing  edge  as  number 
26,  and  the  fictional  exit  station  as  number  30.  The  seven 
stations  prior  to  the  fictional  entrance  station  are  used  to 
generate  the  velocity  profile  expected  at  the  leading  edge.  If 
an  amount  of  prerotation  Is  expected,  the  proper  turning  of  the 
fluid  Is  Imparted  by  a fictional  stator  located  from  stations 
3 to  6.  For  the  specified  amount  of  prerotation,  which  in  this 
case  is  assumed  uniform  across  the  leading  edge,  the  necessary 
tangential  velocity  distribution  at  the  stator  exit  is  calculated. 
Differences  in  radial  location  between  the  stator  exit  and  the 
rotor  leading  edge  are  accounted  for  by  applying  the  equation 
for  free  vortex  flow  (r  Vg  = constant).  The  angular  momentum 
Is  thereby  conserved  since  no  forces  are  acting  on  the  fluid  between 
these  two  stations. 

Within  each  of  the  various  segments  of  the  pump  illustrated 
In  Figure  4,  the  streamlines  are  further  subdivided  into  equal 
increments.  The  impeller  blading  therefore  consists  of  11  stations, 
10  of  which  are  to  be  loaded,  imparting  the  turning  to  the  flow. 

The  output  from  the  Geometry  Module  is  included  in  Appendix  B. 

Along  with  the  coordinates  for  each  of  the  30  stations  and  the 
3 nominal  streamlines,  a listing  is  provided  of  the  additional 
values  needed  as  inputs  by  the  SCM  Module  which  are  calculated 
by  the  GEOMTRY  program.  These  nondimensional  values  are: 

(1)  an  area  mass  flow  coefficient  CMS,  (2)  the  Vg  distribution 
required  at  the  stator  exit  V0S,  (3)  an  advance  ratio  ADJ, 
and  (4)  the  Vg  distribution  required  at  the  rotor  exit,  VSR. 
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where  r^g£  Is  set  equal  to  one  Inch  and  V^g£  is  equal  to  the 
uniform  axial  velocity  at  inlet  V . An  expression  for  U and 

XjC< 

can  be  found  using  Equation  (4);  and  Equation  (19),  which 

liL 

is  derived  from  the  Euler  pump  equation  [Equation  (7)],  can  be 


rewritten: 


V0R 


(20) 


TE 


<Pn 


The  plot  shown  in  Appendix  B is  the  graphical  representation 
of  the  two-dimensional  axlsymmetric  pump.  The  fictional  entrance 
and  exit  stations  are  included,  but  the  sections  which  include 
the  stator  are  omitted.  All  dimensions  given  are  true;  that 
is,  they  are  not  affected  by  the  size  of  the  plot  indicated  by 


the  scale  ratio. 


i 


( ^ 
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V.  STREAMLINE  CURVATURE  METHOD  MODULE . (SCM) 

At  this  point,  the  problem  of  determining  the  axisymmetric 
flow  field  within  the  impeller  channel  can  be  solved.  This 
is  accomplished  by  using  the  computer  program  described  in 
Reference  6.  Other  than  a few  small  changes  needed  to  output 
data  cards,  the  program  remains  intact  and  is  included  in  the 
design  process  as  the  SCM  Module. 

The  Streamline  Curvature  Method  (SCM)  is  able  to  define  the 
streamlines  for  a meridional  surface  and  calculate  the  radial 
and  convective  accelerations  of  the  fluid  based  on  the  stream- 
line geometry.  The  program  performs  an  iterative  procedure 
which  continually  approximates  the  streamline  locations.  It 
is  considered  converged  when  the  maximum  shift  in  any  one  of 
the  streamlines  between  consecutive  passes  is  reduced  to  an 
acceptable  tolerance,  normally  one-tenth  of  one  percent. 

The  major  consideration  required  of  the  designer  at  this 

point  is  how  to  prescribe  the  streamwise  load  along  the  rotor 

blade  stations.  With  the  leading  edge  unloaded,  there  remain 

ten  stations  which  are  required  to  impart  the  necessary  turning 

to  the  fluid  in  order  to  produce  the  V_  distribution  specified 

y 

at  the  rotor  exit  station.  For  the  example  design,  a uniform 
Increase  in  the  angular  momentum  was  selected.  The  percent  of 
the  total  load  at  each  station  for  the  uniformly  loaded  blades 
is  listed  in  Table  1.  It  was  later  determined  that  a different 
loading  distribution  would  be  more  suitable  for  centrifugal  pump 
blading^. 

In  order  to  enhance  the  resistance  to  cavitation  inception 
blades  are  normally  trailing  edge  loaded.  With  the  rotor 
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pressure  distribution  from  Reference  5,  reproduced  in  Figure  5, 
the  static  pressure  on  the  blade  surface  is  not  required  to  rise 
too  abruptly  over  the  first  few  stations  where  the  cavitation 
resistance  is  at  a minimum.  With  an  integration  procedure 
given  in  Appendix  III  of  Reference  5,  a loading  distribution 
for  the  data  from  Figure  5 was  calculated  and  is  included  in 
Table  1.  It  can  be  seen  that  the  percent  of  the  total  load  is 
much  lower  at  first  and  does  not  equal  that  for  a uniformly 
loaded  blade  until  the  seventh  blade  station  (number  22) . From 
that  point  on,  it  clearly  exceeds  the  uniformly  loaded  blade; 
hence,  its  designation  as  trailing  edge  loaded.  The  selection 
of  the  loading  distribution  will  directly  affect  the  velocity 
and  streamline  locations  calculated  by  the  SCM  program  and 
ultimately  the  shape  of  the  blades. 

A sample  output  from  the  SCM  Module  is  listed  in  Appendix 
C.  For  each  station,  the  nondimensional  meridional  and 
tangential  velocities,  the  static  and  total  pressure  coefficients, 
the  streamline  location,  and  the  local  radius  of  curvature  are 
calculated.  The  31  streamlines  used  in  the  example  design  consist 
of  11  major  streamline  locations  with  two  pseudo-streamlines 
Included  between  each  for  a more  accurate  solution.  As  many  as 
51  streamlines  can  be  used,  but  as  an  example,  only  7 streamlines 
were  considered  for  use  in  the  succeeding  modules.  These  are 
identified  as  the  1st,  6th,  11th,  16th,  21st,  26th,  and  31st 
streamline  locations.  For  the  three-dimensional  design  of  the 
blading,  only  stations  8 through  26  are  needed,  the  required 
data  being  obtained  in  the  form  of  punched  cards.  The  data 
In  the  example  were  obtained  after  200  passes  of  the  program 
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where  the  accuracy  had  reached  six-hundredths  of  one  percent 
maximum  shift  in  the  streamline  locations. 
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VI.  STREAMSURFACE  MODULE  (STREAM) 

The  Streamsurface  Module  (STREAM)  begins  the  three-dimensional 

design  of  the  impeller  blading.  Since  it  is  extremely  difficult 

to  handle  three  coordinates  at  any  one  time,  a method  of 

9 

transformation  has  been  employed  . In  each  view  of  the  Impeller, 
one  of  the  coordinates  is  fixed  such  that  a two-dimensional 
representation  of  the  blade  shapes  can  be  constructed.  The 
first  step  in  this  process  is  to  lay  out  the  conformal  grid. 

The  method,  originally  documented  by  Prasil,^  transforms  the 
streamsurface  into  a plane  surface  in  a manner  similar  to  the 
Mercator  Projection  of  the  earth  onto  a map.  The  ordinate  of 
the  conformal  grid  corresponds  to  the  distance  along  the  stream- 
line (the  arc  length  coordinate,  S),  and  the  abscissa  is  the 
blade  wrap  (the  tangential  coordinate,  6).  It  is  the  purpose 
of  the  Streamsurface  Module  to  locate  the  horizontal  grid  lines 
on  the  streamsurfaces  and  determine  the  vertical  position  on 
the  grid  of  each  of  the  blade  stations. 

At  this  point,  it  is  necessary  to  examine  several  figures 
In  order  to  understand  the  three-dimensional  nature  of  the 
problem.  The  conformal  grid  on  which  the  blades  are  to  be 
constructed  is  actually  an  output  from  the  Blade  Shape  Module, 
and  is  therefore  included  in  Appendix  E as  Figure  E-1.  It 
displays  the  S and  0 coordinates  for  each  streamline.  The 
output  of  the  Streamsurface  Module  is  the  plot  shown  in 
Figure  D-1.  It  is  a profile  view  of  the  pump  which  displays 
the  R and  Z coordinates  for  each  streamline.  If  a point  can 
be  determined  by  its  location  on  a given  streamline,  then  the 


radial  and  axial  coordinates  can  be  found. 
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As  an  aid  In  visualizing  the  grid  construction.  Figure  6 
illustrates  the  coordinate  system  for  two  radial  plane  sections 
which  intersect  the  streamsurface.  The  conformal  grid  is 
constructed  for  a given  increment  in  the  tangential  direction 
A6  which  corresponds  to  the  Interval  at  which  radial  planes 
extending  from  the  axis  of  rotation  intersect  the  streamsurface. 
This  increment  is  assigned  a value  of  arc  length  Ax  which 
establishes  the  size  of  the  square  grid  pattern.  A new  reference 
radius  r'  , is  defined  such  that 


(A0)  = Ax 


(21) 


The  law  of  conformal  mapping  then  applies  to  locations  at  radii 
different  than  the  reference  radius.  Any  increment  A5.  on  the 
conformal  grid  that  corresponds  to  a radius  different  than  the 
reference  radius  will  be  proportional  to  the  true  increment 
size  by  the  inverse  of  the  radius  ratio: 


^^CONF  r ^^TRUE 


(22) 


Thus,  the  blade  thickness  at  a radius  less  than  the  reference 
radius  will  appear  larger  on  the  conformal  grid. 

With  a knowledge  of  how  the  conformal  grid  is  to  be 
constructed,  it  is  now  possible  to  establish  the  grid  points 
on  the  meridional  streamsurface.  In  Figure  6 the  first  "square" 
is  identified  by  locating  point  1 on  the  streamline  which 
generates  the  streamsurface.  The  center  of  this  "square", 
point  C,  is  located  on  the  radial  plane  which  corresponds  to 
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an  angle  of  1/2  A0.  The  distance  along  the  STREAMLINES  is 
divided  such  that  the  arc  length  from  point  0 to  the  center 
equals  the  arc  length  from  the  center  to  point  1.  This  section 
approximates  a square  since  the  distance  from  its  center  to 
each  of  Its  sides  Is  equal, 

1/2  r A0  ■ As 

Since  each  streamline  was  defined  by  the  SCM  program,  points 
can  be  established  for  each  horizontal  grid  line,  lines  1-6 
Figure  E~l,  by  locating  the  centers  of  these  "squares"  such 
that  Equation  (23)  holds.  The  result  here  is  that  "squares" 
are  Identified  which  get  smaller  towards  the  leading  edge, 
since  the  arc  length  for  a given  A6  decreases  as  the  radius 
decreases.  (The  left  side  of  Figure  D-1  is  a plot  of  the  arc 
length  as  a function  of  the  radius  for  a given  A0.  It  Is  used 
to  determine  the  location  of  the  grid  points  and  Is  included  here 
only  as  an  aid  In  construction.)  The  locations  of  grid  points 
on  each  streamline  are  Indicated  by  circles  on  Figure  D-1. 

The  locations  on  the  conformal  map  of  each  of  the  stations 
within  the  pump  are  found  by  calculating  the  percentage  of  arc 
length  along  the  streamlines  between  horizontal  grid  lines. 

These  grid  lines  on  the  profile  view  of  the  impeller  are 
represented  by  the  dashed  lines  in  Figure  D-1;  the  crosses 
indicate  the  locations  of  the  stations. 

Included  in  Appendix  D is  the  x-y  coordinate  location  of 
the  intersection  of  several  of  the  streamlines  with  the  grid 


(23) 
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lines  (the  streamsurface  sections).  These  coordinates  correspond 
to  the  axial  depth  Z and  the  radius  R,  respectively.  In  a 
separate  listing,  the  location  of  each  blade  station  Is  given 
In  terms  of  Its  vertical  position  on  the  conformal  grid,  a 
value  of  zero  corresponding  to  the  top  of  the  grid  or  the 
exit  station  at  the  tip  diameter.  These  data  are  obtained  In 
the  form  of  punched  cards  and  used  as  Input  for  the  next  module. 
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Vll.  BLADE  SHAPE  MODULE  (BLADE) 

Once  the  conformal  grid  Is  constructed,  it  is  then  possible 
to  lay  out  the  exact  form  of  the  blade  shape  for  each  of  the 
streamlines.  The  purpose  of  the  Blade  Shape  Module  (BLADE)  is 
to  construct  these  blade  shapes  on  the  conformal  grid,  thereby 
reducing  the  problem  to  a geometrical  transfer  of  coordinates 
from  one  view  of  the  impeller  to  another.  The  first  step  in 
this  process  is  to  determine  the  location  of  the  blade  centerline. 

Before  a decision  is  made  as  to  how  the  blade  centerline 
should  be  placed  on  the  conformal  grid,  the  location  of  the 
predicted  mean  streamline  can  be  determined  by  the  calculation 
of  the  blade  angles.  For  each  station  within  the  impeller,  the 
meridional  and  tangential  components  of  the  fluid  velocity  are 
known  from  the  SCM  Module.  Referring  back  to  the  DESMOD  output 
in  Appendix  A,  the  blade  angle  can  be  determined  from  the  velocity 
diagram. 


tan  6*  = . (24) 

U-  Vg 

The  blade  velocity  U at  each  blade  station  is  calculated  using 
Equation  (11);  the  radial  location  being  known  from  the  SCM 
output. 


U 


(25) 


The  meridional  and  tangential  velocities,  however,  must  be  modified 
to  account  for  real  fluid  effects. 
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Due  to  blade  blockage  and  boundary  layer  growth,  the  bulk 

meridional  velocity  will  be  accelerated.  This  can  be  accounted 

for  by  multiplying  V by  a blockage  factor,  1 + PBL,  which  at 

m 

this  point,  must  be  estimated.  The  example  program  was  run  with 

a uniform  amount  of  blockage  at  each  station  except  for  the 

leading  edge  where  the  blockage  was  assumed  negligible.  For 

the  leading  stations,  where  the  boundary  layer  blockage  is 

minimal,  the  blade  blockage  will  be  significant  due  to  the 

reduced  blade  spacing  at  the  smaller  radii.  However,  for  the 

trailing  stations  where  the  blade  blockage  is  reduced  due  to 

greater  blade  spacing  at  the  larger  radii,  the  boundary  layer 

blockage  becomes  significant.  It  can  readily  be  estimated 

that  these  two  effects  combine  to  yield  a blockage  factor 

2 

which  is  nearly  constant  the  entire  length  of  the  blade. 

Actually,  once  the  blade  shapes  are  constructed,  the  true  blade 

blockage  and  boundary  layer  growth  can  be  calculated.  At  that 

point  the  designer  may  want  to  modify  his  blockage  estimate 

and  run  this  module  over  again. 

The  value  for  Vg  in  Equation  (24)  will  be  different  from 

the  fluid  tangential  velocity  V-  due  to  real  fluid  slip  as 

o 

mentioned  before.  Therefore,  it  is  necessary  to  approximate 
the  deviation  of  the  fluid  angle  from  the  blade  angle  towards 
the  blade  trailing  edge.  A slip  distribution  must  then  be 
estimated  in  a manner  similar  to  that  for  the  blockage 
distribution.  The  program  currently  increases  the  percentage 
slip  from  zero  at  the  leading  edge  to  the  value  specified  at 
the  trailing  edge.  Another  possible  distribution  for  the 
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slip  would  be  to  Increase  Vq  linearly  from  the  fourth  or  fifth 

3 9 

blade  station  to  the  value  specified  at  the  trailing  edge.  ' 

A comparison  of  these  two  methods  is  shown  in  Figure  7. 

Perhaps  the  linear  variation  for  the  slip  is  a bit  more  realistic; 
but,  in  either  case,  the  real  fluid  phenomena  or  the  extent  of 
its  effect  within  the  blade  passage  are  unknown.  However,  at 
the  exit  station  a good  estimate  can  be  made  for  the  slip  factor 
as  defined  by  Equation  (15).  The  methods  of  Stodola  and  Busemann 
as  outlined  in  Reference  8 have  proven  highly  successful  in 
determining  this  factor  at  the  trailing  edge. 

9 

Wisllcenus  has  suggested  that  a value  for  the  estimated 
slip  at  the  blade  trailing  edge  should  be  within  five  percent 
when  assuming  a value  for  1 - PSL  = (Vg/Vg  ) equal  to  0.8.  However, 

It 

he  defines  Vg  as  the  tangential  component  of  the  velocity 
based  on  the  relative  flow,  cross-sectional  area  of  the  blade 

g 

channel  at  the  trailing  edge.  Thus,  the  angle  at  which  the 
fluid  exits  the  blade  channel  is  an  average  between  the  specified 
blade  angle  at  the  trailing  edge  and  the  larger  angle  of  the 
adjacent  blade  at  the  exit  cross-section.  This  difference  is 
Illustrated  in  Figure  8.  Since  this  cross-section  cannot  be 
determined  exactly  until  the  blade  shape  is  completely  designed, 
a value  must  be  estimated  for  the  velocity  ratio  at  this  point. 

The  effect  of  this  approximation  is  to  increase  the  required 
blade  angle  at  exit  by  a few  degrees;  hence,  a value  for  the 
slip  ratio  a few  percentage  points  below  0.8  should  be  used 

^This  procedure  is  based  on  the  Stodola  correction  and  ways  to 
estimate  this  cross-sectional  area,  and  thus  the  velocities, 
are  available. 
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Instead,  (As  can  be  seen  from  Figure  A-1,  smaller  values  of 

A 

Xg  will  result  in  larger  values  of  ^2  •)  For  the  example  design, 
a value  of  0,78  was  used  to  calculate  the  slip  distribution. 

The  final  equation  used  to  calculate  the  blade  angle  at 
each  station  is 


0 * tan 


-1 


(1  + PBL)  V 


m 


U - 


'0 


(1  - PSL) 


(26) 


The  blade  angles  calculated  in  this  manner  will  most  often 
exhibit  an  "S"  shape  characteristic,  typical  of  centrifugal  pump 
blading.  The  blade  angles  tend  to  increase  near  the  leading 
edge  to  compensate  for  the  initial  blade  blockage  of  the  flow, 
and  then  decrease  slowly  to  the  angle  required  at  the  trailing 
edge.  As  a precaution  against  irregularities  in  the  data,  the 
blade  angles  are  smoothed  in  such  a way  that  distributions 
exceeding  cubic  in  nature  are  eliminated.  That  is,  the 
distribution  of  blade  angles  from  leading  edge  to  trailing  edge 
may  only  change  slope  twice. 

The  blade  angles  at  the  trailing  edge  can  remain  as 
calculated,  or,  if  desired,  can  be  averaged  over  all  of  the 
streamlines.  For  a pump  design  with  strictly  radial  exit,  it 
seems  reasonable  to  apply  a constant  blade  angle  at  the  rotor 


The  blade  angles  at  the  leading  edge  should  be  modified 
somewhat  to  account  for  Incidence  of  the  flow.  If  the  blade 
angle  is  in  line  with  the  entering  fluid  angle,  the  blade 
is  said  to  be  at  zero  incidence.  There  may  be  some  justification. 


w 
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however,  to  design  the  blades  to  operate  with  a small  amount  of 
incidence. 

There  appears  to  be  some  Justification  for  designing  the 

blades  with  a small  amount  of  negative  incidence  as  measured 

2 

from  the  blade  centerline.  (See  Figure  9.)  The  effort  here 

is  to  aid  the  design  in  terms  of  cavitation.  Cavitation  will 

occur,  and  a cavity  will  form  on  the  suction  surface  of  the 

blade  when  there  is  an  excessive  amount  of  positive  incidence. 

Even  moderate  values  of  positive  incidence  can  result  in 

excessive  accelerations,  thereby  increasing  the  chance  of 

cavitation  inception.  To  assist  the  flow  in  being  turned  by 

the  blade  over  the  suction  surface,  and  to  guard  against  the 

effects  possible  with  positive  incidence,  it  seems  reasonable 

to  apply  a small  amount  of  negative  incidence.  However,  in 

no  way  should  the  negative  incidence  cause  the  fluid  streamline 

to  extend  beyond  the  pressure  surface.  This  can  result  in 

excessive  drag  on  the  blades  and  may  lead  to  separation,  and/or 

cavitation,  at  the  leading  edge  on  the  pressure  surface. 

A similar  way  of  accomplishing  this  same  effect  is  to  align 

the  pressure  surface  of  the  blade  with  the  mean  streamline  as 

9 

calculated  by  the  blade  angles.  This  in  effect  guarantees 
that  the  flow  does  not  go  beyond  the  pressure  surface  and  still 
helps  the  flow  to  be  turned  by  the  blade.  Further  investigation 
Into  this  matter  may  lead  to  a more  satisfactory  method  to 
handle  the  Incidence  problem  at  the  leading  edge. 

For  the  example  design,  the  centerline  of  the  blade  was 
selected  to  align  with  the  predicted  mean  streamline;  that  is. 


I 

I 


% 
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no  further  modifications  to  the  blade  centerline  were  made.  It 
Is  possible,  though,  to  offset  the  blade  centerline  or  camber 
line  from  the  mean  streamline  as  is  typical  of  blade  design  procedures 
for  axial  machines.^  However,  since  there  is  no  precedent  as 
to  how  much  offset  should  be  incorporated,  it  was  thought 
sufficient  at  this  point  to  allow  the  blade  camber  line  and 
the  predicted  mean  streamline  to  coincide. 

The  effects  of  blockage,  slip,  and  incidence  were  all 
included  in  the  calculation  required  to  place  the  blade  center- 
line  on  the  conformal  grid.  Since  fluid  and  blade  angles 
appear  true  on  this  type  of  mapping  (conformal) , it  is  then 
possible  to  calculate  the  blade  wrap  at  each  station.  To  aid 
in  understanding  this  procedure  a sample  blade  shape  (the  hub 
streamline)  is  constructed  in  Figure  10,  identical  to  that 
appearing  in  Figure  E-1.  The  vertical  location  of  each  station 
from  the  STREAM  output  and  the  calculated  blade  angles  were 
used  to  calculate  the  horizontal  location  of  each  station  on 
the  conformal  grid.  Except  for  stations  2 and  10,  each  segment 
of  the  blade  centerline  was  calculated  with  the  station  location 
as  its  midpoint.  The  first  and  last  segments  (from  stations  1 
to  2,  and  10  to  11)  were  drawn  between  stations  so  that  they 
would  align  with  the  leading  and  trailing  edge  blade  angles. 

Once  the  trailing  edge  is  fixed  on  the  conformal  map  (any 
convenient  location) , the  blade  center  or  camber  line  can  be 
plotted. 

At  this  point,  a thickness  distribution  is  selected  and  the 


pressure  and  suction  surface  plotted.  A constant  blade  thickness 
is  normally  used,  which  in  many  cases  may  be  limited  by 
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foundry  considerations  for  smaller-sized  pumps.  The  example 
design  employed  a constant  blade  thickness  equal  to  the 
trailing  edge  thickness.  Beyond  a point  corresponding  to  a 
specified  amount  of  the  blade  chord  (for  example,  fifteen  percent), 
the  thickness  was  allowed  to  taper  linearly  until  it  reached 
the  leading  edge  thickness.  It  must  be  remembered,  however,  that 
the  thickness  to  be  added  to  the  blade  centerline  must  be  the 
conformal  thickness. 

The  true  thickness  to  be  added  will  thus  be  distorted  due 
to  the  mapping.  A plot  of  the  thickness  distribution  as  a 
function  of  the  radius  in  included  in  Figure  E-1  just  above  the 
grid  lines  and  displaced  to  the  right  such  that  zero  radius 
corresponds  to  vertical  grid  line  D.  The  trailing  edge  thickness 
is  true  at  the  reference  radius,  but  decreases  for  larger  radii 
and  Increases  for  smaller  radii  as  required  by  the  law  of 
conformal  mapping.  The  radial  location  of  each  streamline  will 
then  determine  the  thickness  to  be  added  to  each  blade  station 
in  the  conformal  representation  of  the  blade  shape. 

As  an  aid  in  determining  the  radial  location  of  any  point 
on  the  blade,  a mapping  of  the  radius  versus  the  distance 
along  the  streamline  is  plotted  on  the  grid  for  each  streamline. 
This  mapping  is  constructed  using  the  known  radial  locations 
and  vertical  grid  coordinates  for  each  of  the  blade  stations 
and  the  horizontal  grid  lines.  Once  a point  on  the  blade 
centerline  is  identified,  the  radius  can  be  determined  by 
moving  horizontally  on  the  grid.  Figure  E-1,  intersecting 
the  radius  mapping  for  that  particular  streamline.  The 
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thickness  at  that  point  can  then  be  determined  by  moving 
vertically  on  the  grid,  intersecting  the  thickness  mapping, 
j If  this  procedure  is  followed  for  each  station  of  each 

i streamline,  the  three  surfaces  of  the  blade  shape  (suction, 

i pressure,  and  centerline)  can  be  plotted  on  the  conformal 

‘ grid.  The  blade  shape  is  completed  by  constructing  a circle 

I 

for  the  leading  edge  with  a radius  equal  to  half  of  the 
leading  edge  thickness.  It  is  also  possible  to  use  other 
leading  edge  shapes.^ 

' At  this  point,  the  design  of  the  blading  is  complete.  It 
is  now  only  a problem  of  transferring  these  conformal  blade 
shapes  just  constructed  to  the  proper  coordinates  required  by 
each  view  of  the  impeller.  The  objective  here  is  to  construct 
board  sections'*^  on  the  axial  or  plan  view  of  the  impeller.  A 
plot  of  this  nature,  similar  to  the  vay  in  which  hull  shapes 
are  specified  by  the  naval  architect,  allows  the  blading  to  be 
fabricated  by  the  pattern  maker.  This  fabrication  process 
described  in  Reference  1 makes  use  of  boards,  cut  to  the 
shape  of  the  blade  and  of  a uniform  thickness,  which  correspond 
to  sections  passed  perpendicular  to  the  axis  of  rotation  and 
Intersecting  the  blade  surface.  These  can  be  glued  together 
and  smoothed  to  obtain  the  exact  surface  specified  by  the 
design.  Metal  blades  could  then  be  cast  if  the  wooden  blade 
is  used  to  make  a mold.  It  is  also  possible  to  numerically 
mill  the  blades  out  of  metal  since  the  milling  machine  is 
capable  of  cutting  the  board  sections  directly. 

"^^These  will  be  constructed  by  the  Board  Section  Module  and  are 
Illustrated  in  Figure  G-1. 
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In  order  to  construct  the  board  sections  in  the  axial  view 
of  the  impeller,  the  radial  and  tangential  (angular)  coordinates 
(R  and  0)  are  needed  for  each  board  section  location  or  axial 
depth  (Z  coordinate).  From  the  conformal  map,  it  is  possible 
to  obtain  both  the  R and  0 coordinates,  but  no  information  can 
be  found  as  to  the  Z coordinate.  It  is  therefore  necessary  to 
transfer  the  blade  shapes  to  the  profile  or  side  view  of  the  impeller, 
using  the  R coordinate  and  the  streamline  location,  for  a given 
radial  plane  or  0 coordinate.  In  this  way,  sections  of  known 
tangential  location  can  be  Identified  and  the  R and  Z coordinates 
obtained.  These  sections,  derived  from  passing  planes  parallel 
to  the  axis  of  rotation,  are  labeled  radial  sections  and  will 
be  constructed  by  the  next  module  RADIAL.  The  Blade  Shape  Module 
is  therefore  required  to  output  all  of  the  coordinate  locations 
needed  as  input  by  the  Radial  Section  Module. 

Each  vertical  line  (lines  A through  G)  on  the  conformal 
grid.  Figure  E-1,  corresponds  to  a different  tangential  location. 
Typically,  several  lines  (which  are  not  shown  in  Figure  E-1) 
will  be  passed  between  vertical  grid  lines,  thereby  determining 
additional  radial  sections  for  greater  accuracy.  The  radial 
coordinate  for  each  tangential  location  can  then  be  determined, 
and  these  coordinates  for  a sample  streamline  are  included  in 
Appendix  E as  Table  E-1.  The  trailing  edge  (which  is  not  a true 
radial  section)  is  identified  by  the  entry  in  the  last  column 
as  radial  section  number  1,  and  its  centerline  is  usually 
positioned  such  that  grid  line  A is  radial  section  number  2. 

A problem  arises,  however,  in  that  the  hub  streamline  blade 
shape  will  require  more  radial  sections  than  the  tip  streamline 


d 
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blade  shape.  A way  must  therefore  be  found  to  handle  this 
leading  edge  problem. 

Because  the  leading  edge  Is  not  of  Itself  a streamline, 
coordinate  locations  along  it  are  only  known  at  discreet 
points  (i.e.,  where  the  leading  edge  Is  Intersected  by  the 

j 

streamlines).  Obviously  for  an  infinite  number  of  streamlines, 
the  leading  edge  would  be  completely  characterized  by  the 
streamline  stations.  However,  with  only  a finite  number  of 
these  locations  known,  the  leading  edge  has  to  be  mapped 
separately  in  order  to  transfer  it  from  one  view  to  the  next. 

For  reasons  which  will  become  clear  in  the  discussion  of 
the  next  module,  the  leading  edge  was  characterized  by 
identifying  nine  distinct  points  on  each  streamline.  Thus, 
the  leading  edge  now  consists  of  nine  lines,  some  of  which  nay 
not  be  Independent  in  any  one  view.  Radial  sections  which 
would  intersect  streamlines  beyond  the  leading  edge  in  Figure 
E-1,  can  now  be  seen  to  intersect  these  leading  edge  lines.  By 
identifying  the  leading  edge  lines  on  the  profile  view  of  the 
impeller,  incomplete  radial  sections  which  do  not  Intersect  all 
of  the  streamline  blade  shapes  can  be  completed  if  the  radial 
coordinate  is  obtained  for  the  intersection  of  the  radial  section 
with  each  leading  edge  line. 

In  order  to  find  the  radial  coordinate  for  the  leading 
edge  points,  the  streamline  radius  mapping  must  be  used. 

Leading  edge  lines  can  be  located  on  this  auxiliary  mapping 
such  that  the  R coordinate  for  any  point  on  the  leading  edge 
line  can  be  found.  The  R and  0 coordinates  are  listed  in 


Appendix  E for  one  of  the  nine  leading  edge  points  determined 
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by  the  Incomplete  radial  sections.  The  coordinates  for  each 
of  the  leading  edge  points,  Identified  by  crosses  on  the 
streamline  blade  shapes.  Figure  E-1,  will  be  used  to  map  the 
leading  edge  lines  In  both  the  profile  and  the  plan  view  of  the 
Impeller. 


It  is  now  possible  to  evaluate  the  effects  of  stall  or 
separation  on  the  newly  designed  blade  shapes.  The  parameter 


which  will  be  used  for  this  purpose  is  the  lift  coefficient. 
For  a known  value  of  the  retardation  of  the  relative  flow,  a 
limiting  value  of  the  lift  coefficient  can  be  found  from 
Reference  8.  The  suggested  maximum  lift  coefficient  Is 
illustrated  in  Figure  11. 

The  lift  coefficient  Is  calculated  based  on  the  average 
relative  velocity  from  Inlet  to  exit. 


where 


2 ir  r 


W = 1/2  + W^) 


W = + (U  - 


If  a value  Is  selected  for  the  number  of  blades  z,  the  lift 
coefficient  for  each  streamline  can  be  found.  The  value  to  be 
used  for  the  chord  length,  however,  is  somewhat  deceptive. 
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(Note  that  the  use  of  the  conformal  chprd  length  would  be 
Incorrect.)  The  true  chord  length  is  the  distance  covered  by 
the  fluid  as  it  flows  along  the  streamline.  For  Francis-type 
Impellers,  not  only  does  the  streamline  distance  increase  in 
the  tangential  and  axial  directions,  but  it  also  Increases  in 
the  radial  direction.  From  the  axisymmetric  flow  field,  the 
distance  along  the  streamline  is  known  (i.e.,  changes  in  R and 
Z),  but  the  tangential  distance  covered  must  be  approximated. 
This  was  done  by  taking  very  small  increments  in  length  along 
the  streamline,  assuming  that  the  chord  length  could  be 
approximated  for  each  increment  as  a straight  line.  Figure  12 
Illustrates  the  method  used  to  calculate  the  true  chord  length. 
In  this  calculation,  segments  of  chord  length  were  equated  to 
straight  line  segments  which  connect  the  two  sides  of  a right 
triangle  — the  distance  along  the  streamline  As,  and  the  blade 
wrap  R A9. 

The  lift  coefficient  was  found  for  each  streamline  using 
several  different  numbers  of  blades.  The  values  for  ten  blades 
are  outputed  by  the  Blade  Shape  Module  and  are  listed  in 
Appendix  E.  It  can  be  seen  that  the  critical  streamline  is 
the  tip  streamline  (number  1)  where  the  chord  length  is  usually 
a minimum,  thus  producing  a maximum  lift  coefficient  at  the 
smallest  value  for  the  retardation  of  the  relative  flow. 

Based  on  these  values  the  designer  can  select  the  number 


JL 


of  blades  to  be  used  in  the  actual  pump  by  making  use  of 
Figure  11.  Care  should  be  taken  in  that  values  for  the  lift 
coefficient  close  to  the  suggested  limit  may  not  be  as 
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conservative  as  the  designer  might  wish  them  to  be.  If  the 
selection  of  the  number  of  blades  at  this  point  Is  not 
satisfactory,  he  could  fix  the  number  of  blades  and  redesign 
the  Impeller  profile  to  obtain  better  values  for  the  lift 
coefficient.  The  lift  coefficient  will  be  reduced  for  blades 
of  longer  chord  length  by  locating  the  leading  edge  closer  to, 
or  locating  the  exit  station  farther  (axially)  from  the 
entrance  station.  The  value  for  the  retardation  of  the 
relative  flow  will  be  Increased  for  designs  of  larger  tip 
diameter  and/or  smaller  exit  station  width  (assuming  the 
parameters  at  the  entrance  station  are  fixed). 
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VIII.  RADIAL  SECTION  MODULE  (RADIAL) 

In  order  to  transfer  blade  shapes  to  the  profile  view  of 
the  Impeller,  the  streamline  number  and  the  radial  coordinate 
must  be  known.  This  Information  was  obtained  by  the  previous 
module  and  enables  the  program  to  construct  radial  sections 
(which  are  determined  by  a given  radial  plane)  In  this  view. 
These  sections  are  plotted  Identifying  the  suction,  pressure, 
and  centerline  surfaces.  It  Is  the  purpose  of  the  Radial 
Section  Module  (RADIAL)  to  perform  this  transfer  and  determine 
the  board  section  coordinates  required  as  Input  by  the  next 
module. 

The  BLADE  output  data  can  be  plotted  directly  and  Is  shown 
In  Figure  F-1  of  Appendix  F.  Radial  sections  2 through  26  are 
seen  to  span  the  blade  channel  from  the  hub  streamline  to  the 
tip  streamline.  The  radial  sections  which  correspond  to  a 
vertical  grid  line  on  the  conformal  map  (Figure  E-1)  are 
Identified  by  letters  A through  G.  Also  shown  in  Figure  F-1 
are  the  vertical  board  sections  numbered  1 to  19.  Up  to  the 
20th  radial  section,  which  corresponds  to  a tangential  location 
on  the  conformal  grid  midway  between  grid  lines  E and  F,  the 
radial  sections  are  complete  in  that  the  surfaces  Intersect 
all  seven  streamlines.  The  pressure  surface  of  these  radial 
sections  Is  the  upper  solid  line  and  the  suction  surface 
Is  the  lower  solid  line.  [From  the  conformal  grid  (Figure  E-1), 
it  can  be  seen  that  a vertical  line  (or  radial  section) 
intersects  the  suction  surface  (the  lower  one)  at  a point 
closer  to  the  leading  edge  of  the  blade  shape,  which,  from 
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the  radius  mapping,  corresponds  to  a smaller  radius  than  that 
of  the  centerline.  The  opposite  is  true  for  the  pressure 
surface.]  A problem  does  not  arise  until  the  surfaces  of  the 
radial  sections  fail  to  intersect  all  of  the  streamlines,  thereby 
intersecting  the  leading  edge. 

By  using  the  nine  distinct  leading  edge  points.  Incomplete 
radial  sections  which  Intersect  the  leading  edge  can  be 
drawn  with  the  required  accuracy.  The  only  way  a point  can 
be  identified  on  this  profile  view  is  by  its  streamline  number 
and  radial  coordinate  (R).  Since  the  Incomplete  radial 
sections  terminate  at  locations  between  streamlines,  they 
cannot  be  identified  directly. 

The  nine  leading  edge  points  are  defined  from  the  conformal 
map  and  illustrated  in  Figure  13.  Points  are  identified  on 
each  of  the  suction,  pressure,  and  centerline  surfaces  at  four 
tangential  (or  radial  plane)  locations.  The  "c"  location  is 
calculated  to  be  the  point  of  maximum  tangential  position  and 
as  such  can  only  intersect  the  pressure  surface.  The  "b" 
location  Intersects  the  centerline  at  its  maximum  tangential 
position  and  therefore  does  not  intersect  the  suction  surface. 

The  "a"  location  corresponds  to  the  center  of  the  leading 
edge  circle  and  also  Identifies  the  point  of  minimum  radius  on 
the  suction  surface.  The  "z"  location  is  determined  such  that 
the  "z"  pressure  surface  point  is  at  the  same  radius  as  the 
"b"  centerline  point.  It  should  be  clear  that  only  in  the 
conformal  map,  where  the  blade  is  constructed,  is  the  leading 
edge  actually  defined  by  a circle.  This  leading  edge  shape 


-45- 


5 January  1978 
MCB:jep 


will  change  from  view  to  view,  depending  on  the  coordinates 
of  the  leading  edge  points  for  the  various  types  of  sections. 

The  leading  edge  points  identified  on  the  conformal  grid 
can  be  mapped  onto  the  profile  view  of  the  impeller.  These 
leading  edge  lines,  plotted  on  Figure  F-1,  represent  the  area 
in  which  the  leading  edge  shapes  of  the  Incomplete  radial 
sections  are  to  be  plotted.  (Because  of  the  reduced  scale  in 
Figure  F-1,  the  leading  edge  lines  tend  to  obscure  the  shape 
of  the  sections  within  the  leading  edge  area.)  A typical 
leading  edge  configuration  is  shown  plotted  in  Figure  14.  The 
grid  in  this  figure  consists  of  blade  surfaces  intersected  by 
streamlines;  however,  these  streamlines  are  not  necessarily 
those  which  are  already  known.  Upon  examination  of  the 
conformal  mapping.  Figure  E-1,  one  will  find  that  a radial 
section  (or  vertical  line)  Intersects  the  streamline  blade 
shapes  in  normally  three  places,  one  for  each  surface.  Although 
these  points  each  correspond  to  a different  radius  (as  determined 
by  the  radius  mapping),  they  are  all  on  the  same  streamline. 

If  these  three  leading  edge  points  for  each  streamline  are  connected 
OQ  the  conformal  grid  to  form  three  leading  edge  lines  (samples  of 
which  are  shown  in  Figure  13),  a radial  section  which  Intersects 
these  lines  will  necessarily  define  a different  streamline  by 
the  points  of  intersection.  Thus,  vertical  sections  may  be 
taken  through  the  leading  edge  which  will  form  streamlines  parallel 
to  the  known  streamlines.  Since  the  R coordinate  for  each  leading 
edge  point  was  outputed  by  the  Blade  Shape  Module  for  each 
incomplete  radial  section,  the  leading  edge  shapes  of  these 
radial  sections  can  be  plotted. 
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At  this  point,  vertical  sections  can  be  passed  through 
the  radial  sections  plotted  on  the  profile  view  of  the  impeller 
In  Figure  E-1.  These  sections  are  located  using  a given 
Increment  In  the  axial  direction  and  identified  as  board  sections, 
where  the  blade  surfaces  are  Intersected  by  planes  perpendicular 
to  the  axis  of  rotation.  The  radial  and  axial  coordinates  (R 
and  Z)  can  be  found  for  each  board  section  that  Intersects  a 
radial  section  surface.  These  coordinates  are  calculated  by 
the  program  for  each  radial  section;  a sample  output  is  included 
In 'Appendix  F.  Coordinates  are  listed  for  board  sections  that 
intersect  each  particular  radial  section  surface,  beginning 
at  the  hub  streamline  and  continuing  until  the  tip  streamline 
is  reached.  [For  the  board  section  plot,  the  centerline  surface 
Is  not  required.]  Board  sections  which  intersect  the  leading 
edge  lines  will  be  Incomplete  in  the  axial  or  plan  view  of  the 
Impeller.  However,  because  of  the  leading  edge  mapping,  the 
coordinates  can  be  obtained  for  each  board  section  that  inter- 
sects these  leading  edge  lines.  A sample  list  of  these 
coordinates  is  also  included  In  Appendix  F.  Once  the  leading 
edge  points  are  mapped  on  the  plan  view  of  the  Impeller,  the 
Incomplete  board  sections  can  be  constructed  using  these  leading 
edge  coordinates. 

Because  each  blade  shape  is  aligned  with  the  mean  stream- 
line, some  waviness  in  the  radial  sections  of  Figure  F-1  may 
result.  This,  In  and  of  Itself,  is  not  cause  for  alarm,  but 
there  may  be  some  question  as  to  the  shape  of  these  radial 
sections  where  they  intersect  the  shrouds.  Particularly  at  the 
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outer  shroud  (or  tip  streamline).  It  has  been  pointed  out  * 
that  acute  angles  at  this  point  can  result  in  hydrod3mamically 
unfavorable  conditions  within  the  flow  passage.  The  most 
direct  way  to  alleviate  this  condition  is  to  properly  position 
the  leading  edge  at  the  outer  shroud  as  far  back  as  possible, 
and  the  leading  edge  at  the  inner  shroud  as  close  to  the 
entrance  as  practical.  This  will  result  in  a leading  edge 
which  is  closer  to  a concial  surface  than  a cylindrical 
surface.  To  avoid  excessively  long  chord  lengths  at  the  hub 
streamline,  the  exit  blade  station  should  probably  be 
chamfered  (creating  a blade  design  which  more  closely 

resembles  a mixed  flow  impeller) . 

9 

Wislicenus  has  indicated  that  is  may  be  sufficient  to 

design  only  the  tip  streamline  blade  shape  to  align  with  the 

/ 

mean  streamline.  The  tip  streamline  is  the  most  critical  in 
terms  of  cavitation,  separation,  and  frictional  losses.  He 
suggests  that  the  remaining  streamlines  be  designed  and 
positioned  based  on  primarily  geometrical  considerations. 

This  will  help  in  lessening  the  problem  of  acute  angles  at 
the  shrouds,  but  may  result  in  blade  shapes  that  do  not 
impart  the  proper  turning  to  the  flow.  Further  investigation 
is  needed  in  this  area  to  guide  the  designer. 


I 


-48- 


5 January  1978 
MCB: jep 


Ft 

IX.  BOARD  SECTION  MODULE  (BOARD) 

The  Board  Section  Module  (BOARD)  is  the  last  segment  of  the 
computer  program;  the  only  output  from  this  module  is  the 
plot  of  the  plan  or  axial  view  of  the  Impeller.  At  this  point, 
practically  all  of  the  data  points  required  are  known,  the 
majority  of  them  calculated  in  the  previous  module.  The  Board 
Section  Module  is  therefore  only  required  to  plot  the  R - 6 
projection  of  the  blading  surfaces  and  to  identify  the  board 
sections. 

A short  procedure  which  was  omitted  from  the  last  module 
was  incorporated  into  the  BOARD  Module.  Referring  back  to 
Figure  F-1,  it  can  be  seen  that  board  sections  (the  vertical 
lines)  which  intersect  the  hub  to  tip  streamlines  terminate  at 
these  locations.  These  board  sections  are  not  complete  in  that 
they  do  not  extend  from  the  leading  edge  to  the  trailing  edge. 
Where  they  intersect  either  the  hub  or  tip  streamline,  the 
surfaces  of  these  board  sections  are  determined  by  a single 
point  or  radial  location  in  Figure  F-1.  When  these  end  points 
are  transferred  to  the  plan  view,  they  will  be  seen  as  an  arc, 
equivalent  to  the  width  of  the  blade  at  that  point,  and  located 
at  a constant  radius  from  the  axis  of  rotation.  At  first,  board 
sections  will  be  plotted  which  start  at  the  trailing  edge  and 
terminate  on  the  hub  streamline.  Depending  on  the  shape  of  the 
impeller,  a few  board  sections  might  extend  from  the  trailing 
to  the  leading  edges.  The  remaining  board  sections  will  begin 
at  the  tip  streamline  and  terminate  at  the  leading  edge. 

The  cylindrical  coordinate  grid  for  the  plan  view  of  the 
impeller  is  laid  out  in  Figure  G-1,  identifying  the  radial 
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planes  which  appeared  as  vertical  lines  on  the  conformal  grid 
(Figure  E-1).  The  entrance  station  hub  and  tip  radii  are 
plotted  to  provide  the  designer  with  the  view  that  he  would 
see  if  he  looked  directly  into  the  eye  of  the  impeller. 

Since  the  actual  blading  only  extends  from  the  inner 
shroud  to  the  outer  shroud,  the  hub  and  tip  streamlines  will 
appear  in  this  view  as  the  side  boundaries  within  which  the 
board  sections  are  drawn.  Obviously,  the  leading  and  trailing 
edges  represent  the  other  two  boundaries.  The  required  data 
for  these  two  streamlines  (which  are  not  true  board  sections, 
but  the  locus  of  the  endpoints  of  an  infinite  number  of  board 
sections)  were  obtained  by  the  Blade  Shape  Module  from  the 
conformal  grid.  These  streamlines  can  be  plotted  and  the  arcs, 
where  the  board  sections  intersect  either  the  hub  or  tip 
streamlines,  drawn. 

Board  sections  are  drawn  in  this  view  of  the  impeller 
with  the  suction  surface  in  the  foreground  as  a solid  line 
and  the  pressure  surface  behind  it  as  a dashed  line.  These 
surfaces  can  be  plotted  directly  from  the  data  outputed  by  the 
Padial  Section  Module.  Incomplete  board  sections  which  inter- 
sect the  leading  edge  can  be  drawn  only  up  to  the  point 
where  the  pressure  and  suction  surface  begin  to  converge. 

To  aid  in  constructing  the  leading  edge  shapes  for  each 
of  the  incomplete  board  sections,  a mapping  of  the  leading 
edge  lines  on  the  R - 6 grid  is  included  in  Figure  G-1,  but 
displaced  a nominal  distance  in  the  horizontal  direction.  A 
single,  solid  leading  edge  line  (L.E.  point  BP)  is  drawn  to 
connect  the  hub  and  tip  streamlines.  It  represents  the  extreme 
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surface  of  the  blading  as  seen  by  the  designer  looking  into 

the  eye  of  the  impeller.  The  true  or  nominal  leading  edge 

line  (L.E.  point  BC)  is  also  drawn  in  Figure  G-1  as  a dashed 

line;  but,  due  to  the  reduced  size  of  the  plot,  it  is  not  | 

clearly  visible.  j 

The  data  for  the  Incomplete  board  sections  consists  of  j 

a board  section  number  (or  axial  coordinate  Z)  and  a radial  | 

location  (or  R coordinate)  which  corresponds  to  the  point  of 
intersection  of  the  board  section  with  each  of  the  leading 
edge  lines.  The  0 coordinate  can  then  be  determined  for  each 
of  the  nine  leading  edge  points  by  making  use  of  the  leading 
edge  line  mapping.  In  this  way,  the  board  section  leading 

edge  surfaces  are  constructed  in  much  the  same  manner  as  that  | ■ 

used  for  the  incomplete  radial  sections.  Thus,  the  drawings 


for  the  pump  impeller  design  are  completed. 
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X.  SUMMARY  AND  CONCLUSIONS 

The  fabrication  of  the  blading  can  now  be  visualized  more 
clearly.  In  Figure  G-1,  each  board  section,  from  the  large 
sections  near  the  hub  to  the  small  sections  near  the  tip,  can 
be  cut  from  wooden  boards  of  uniform  thickness  corresponding 
to  the  distance  between  the  vertical  board  sections  in 
Figure  F-1.  These  sections  are  laid  up  on  top  of  one  another 
and  staggered  to  approximate  the  blade  shape.  They  are  then 
glued  and  later  smoothed  until  the  true  shape  which  appears 
in  Figure  G-1  is  obtained.  (Certainly  many  more  board  sections 
than  those  illustrated  could  be  used  to  reduce  the  necessary 
smoothing . ) 

Also,  since  the  exact  location  of  the  blade  surface  is 
known  in  the  plan  view  of  the  impeller,  it  is  now  possible  to 
evaluate  the  slip  correction  approximated  earlier.  The  number 
of  blades  defines  the  spacing  at  the  trailing  edge;  and,  by 
reproducing  this  same  blade  shape  at  the  proper  tangential 
Interval,  the  cross-section  at  the  blade  exit  can  be  derived. 

Either  the  methods  of  Stodola  and  Busemann  or  Wlslicenus,  as 
discussed  in  Section  VII,  could  then  be  applied. 

The  only  design  quantities  from  the  original  one-dimensional 
analysis  which  remain  unknown  at  this  point  are  the  amount  of 
leakage  and  the  hydraulic  efficiency.  The  leakage  cannot  be 
determined  until  the  casing  and  seals  are  designed.  This  would 
require  that  the  inlet  housing  as  well  as  the  exit  housing  (normally 
a volute)  be  designed,  which  is  a separate  part  of  the  overall 
pump  design  process. 
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The  hydraulic  efficiency  is  the  true  measure  of  the 
performance  of  the  pump,  yet  despite  this  detailed  analysis. 

It  must  remain  only  an  estimate.  Impellers  designed  in  this 
manner  have  been  built  and  .tested  with  very  good  results."*" 

The  design  techniques  used  by  the  computer  program  are  based 
on  sound  hydrodynamic  principles.  The  hydraulic  efficiency 
has  always  been  determined  by  experimental  testing  for  any 
type  of  pump  regardless  of  its  design.  However,  because  of 
the  detailed  approach  used  here,  there  may  be  a way  to  approximate 
the  efficiency  in  terms  of  fluid  losses.  Since  the  channels 
within  the  blade  row  are  defined  at  many  points  from  the  leading 
edge  to  the  trailing  edge,  the  pertinent  velocity  and  pressure 
data  can  be  determined.  It  may  then  be  possible  to  calculate 
the  frlctiona] , boundary  layer,  and  secondary  losses  in  this 
manner.  Such  an  efficiency  will  of  course  be  only  an 
estimate  and  subject  to  experimental  verification. 

Once  the  designer  is  satisfied  with  the  pump  at  its  design 
operating  condition,  the  next  step  would  be  to  evaluate  its 
off-design  performance.  Here,  it  is  necessary  to  determine 
the  actual  flow  based  on  the  designed  blade  shapes  and  a 
specified  inlet  condition.  Because  of  the  detailed  design  of 
the  blades,  and  the  thorough  characterization  of  the  flow  through 
the  Impeller,  techniques  might  be  derived  to  solve  this,  the 
direct,  problem  of  turbomachinery. 


Wisllcenus  has  been  responsible  for  the  design  of  many  pumps  over 
the  past  40  years.  Perhaps  the  most  famous  of  which  is  the  units 
Installed  at  the  Aswan  Colorado  State  Dam  which  have  demonstrated 
efficiencies  of  over  90  percent. 
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The  user  of  this  pump  impeller  design  computer  program 
must,  in  fact,  be  the  designer  or  someone  working  directly 
under  his  guiding  hand.  The  program  of  itself  is  not  a complete 
design  method,  but  rather  a tool  to  be  used  by  the  designer. 

There  are  may  points  along  the  way  which  require  engineering 
input  and/or  evaluation,  parts  of  which  may  be  altered  to 
incorporate  or  investigate  other  design  techniques.  In  this 
way,  the  designer  is  able  to  interact  with  the  computer  and 
produce  a detailed  hydrodynamic  design  with  which  can  be 
satisfied,  in  a relatively  short  amount  of  time. 

The  basic  advantage  of  this  method  is  its  ability  to 
handle  the  three-dimensional  design  problem  by  making  use 
of  the  various  two-dimensional  mappings.  The  exact  procedures 
programmed  are  a result  of  the  particular  design  method  selected; 
however,  the  pump  design  would  still  proceed  in  much  the  same 
manner  if  another  method  was  used  instead.  Thus,  the  designer 
has  not  only  a tool  which  can  be  used  to  quickly  and  easily 
design  pump  Impellers,  but  also  the  means  to  evaluate  the  design 
by  observing  the  computer  results  of  the  various  design  changes 
that  could  be  made. 

The  goal  of  the  designer  is  to  be  able  to  produce  a pump 
which  will  meet  the  requirements  specified  by  the  system  of 
which  it  is  a part.  This  can  best  be  accomplished  if  the 
techniques  are  available  to  analyze  each  step  of  the  design 
process  in  an  attempt  to  fully  understand  what  is  taking  place. 

It  is  hoped  that  this  pump  Impeller  design  computer  program 
is  a step  in  that  direction. 
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TABLE  1 


ROTOR  BLADE  LOADING  DISTRIBUTIONS 


Station  Percent  of  the  Total  Load 

Uniformly  Loaded  Trailing  Edge  Loaded 


16 

0 

0.0 

17 

10 

6.3 

18 

20 

15.0 

19 

30 

25.2 

20 

40 

36.5 

21 

50 

48.3 

22 

60 

60.5 

23 

70 

72.7 

24 

80 

84.3 

25 

90 

94.5 

26 

100 

100.0 
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Mixed  Flow  Impeller 


Francis-Type  Impeller 


Unconventional  Axial  Flow  Impeller 


Figure  1 - Types  of  Impellers 
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Figure  3 
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- Relative  Configuration  of  the  Inlet  and  Exit  Stations 
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TANGENTIAL  VELOCITY,  (fpsi 


Figure  7 - Approximations  for  the  Deviation  of  the  Fluid 
Angle  from  the  Blade  Angle 
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Figure  8 - Method  to  Determine  the  Exit  Blade  Angle 
to  Account  for  Fluid  Slip 
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Figure  9 - Incidence  Effects  at  the  Leading  Edge 
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Figure  10  - Detailed  Construction  of  the  Blade  Shapes 
on  the  Conformal  Grid 
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Figure  12  - Determination  of  the  True  Chord  Length 
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Figure  14  - Location  of  the  Nine  Leading  Edge  Points  on  Streamlines 
Determined  by  Radial  Sections 
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APPENDIX  A 


DESIGN  MODULE  OUTPUT 
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Table  A-1  (continued) 


Figure  A-1  - Design  Module  Velocity  Diagrams 
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APPENDIX  B 


GEOMETRY  MODULE  OUTPUT 
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Table  B-1  - Impeller  Station  Coordinates 


Table  C-1  - Sample  SCM  Station  Output 
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APPENDIX  D 


STREAMSURFACE  MODULE  OUTPUT 
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Figure  D-1  - Construction  of  the  Conformal  Grid  on  the  Streamsur faces 
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Table  D-2  - Station  Locations  for  Each  Streamline 
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APPENDIX  E 


BLADE  SHAPE  MODULE  OUTPUT 


Table  E-1  - Radial  Section  Coordinates  for  a Sample  Streamline 
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Table  E-1  (continued) 


Table  E-1  (continued) 
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Figure  F-1  - Radial  Sections  Transferred  to  the 
Profile  View  of  the  Impeller 


Table  F-1  - Board  Section  Coordinates  for  a Sample  Radial  Section 
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BOARD  SECTION  MODULE  OUTPUT 


Figure  G-1  - Board  Sections  Jn  the  Plan  View  of  the  Impeller 
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